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Abstract 

Two aspects of hadronic B decays are investigated. Firstly, the supersym- 
metric enhancement of hadronic b decays b — ► qg — ► qq'q', q G {d, s} by gluino 
penguin processes is studied through their effect on the Wilson coefficients of 
the effective Hamiltonian. The gluino penguin is dominated by the magnetic 
dipole transition which is strongly magnified relative to the electric monopole 
driven standard model gluon penguin by the renormalization-induced QCD 
corrections, resulting in quark decay rates for pure penguin processes which, 
at scales 0(mb), can greatly exceed the standard model rates. The CP asym- 
metries are, however, relatively unaffected. Secondly, hadronization of the 
final state quarks is studied through a simple phase space spectator model. 
We consider two extreme models for color flow during meson formation; one 
in which color flow is ignored and one of color suppression in which low mass 
meson formation occurs only for color singlet quark antiquark pairs. We 
find that processes in which the spectator antiquark q s combines with q' are 
relatively insensitive to the color flow model whereas processes in which q s 
combines with q are very sensitive to color suppression. 

12.15.Ji, 12.60. Jv, 13.25.Hw 
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I. INTRODUCTION 



These are exciting times for B physics. Recently the CLEO collaboration |I| reported 
the first measurements of a number of exclusive charmless hadronic B decays which provide 
conclusive evidence for hadronic penguin processes, and the next decade will see || inten- 
sive investigation of the S-meson system at CESR, the Tevatron, HERA, the SLAC and 
KEK .B-factories, and at the LHC. In particular, measurements of rare flavor-changing B 
decays [[3] will provide windows for the discovery of indirect effects of new physics beyond 
the standard model (SM). For example, data on CP asymmetries should indicate whether 
the Cabbibo-Kobayashi-Maskawa (CKM) CP violation mechanism of the SM is correct or 
whether new sources of CP violation are needed. The most favoured candidate theory for 
this new physics is supersymmetry (SUSY) for which rare B decays offer possible insight be- 
cause 0|| the measured observables involve SM and SUSY processes occurring at the same 
order of perturbation theory. In particular, the b — > s transition provides an opportunity 
to study CP violation from non-standard phases ||. Probing for SUSY in b — > sj decays 
has been examined by many authors and there is also significant current interest in 

the b — > sg penguin decay fll3|-|T6[| for which it has been argued [|T7| that enhancement for 
on-shell gluons is needed from non-SM physics to explain the CLEO measurement of a 
large branching ratio for the inclusive process B — > r]'X s and the rf — g — g gluon anomaly. 

The most predictive of the SUSY models is the (constrained) minimally supersymmetric 
standard model (MSSM) [ [TB"| , |H| based on spontaneously broken N = 1 supergravity with flat 
Kahler metrics pO|l , universal explicit soft-SUSY breaking terms at the scale Mmssm ~ Mu 
and spontaneous breaking of the SU(2) ®U(1) symmetry driven by radiative corrections. 
The unification scale Mu boundary conditions are a universal scalar mass m , universal 
gaugino mass mi/2 and a universal soft SUSY-breaking trilinear scalar coupling A. After 
minimization of the full one-loop Higgs effective potential the MSSM is represented by a 
four- dimensional parameter space {m , my2, A, tan/3}, where tan/3 = vi/v\ is the ratio of 
the vacuum expectation values of the two Higgs fields, together with the sign of the coupling 
\x between the two Higgs fields. 

In an earlier paper E1J we examined the SM and MSSM predictions at the weak scale 



quark level for the penguin mediated decays b — > qq'q', q G {d, s}. For the SM we presented 
the relative contributions of the internal u, c and t quarks to the gluon penguin and stressed 
the invalidity of the widespread assumption that the process is dominated by the t quark 
22fl . We also considered the relative magnitudes of the various form factors and the role of 



the strong and weak phases p3| , |2~i|| and found, for example, that the CP violating phases 
for the b — > dg and b — > sg electric form factors, which dominate the decay amplitude, 
have no simple relationship with any angle of the unitarity triangle. For the MSSM we 
explored the allowed regions of the parameter space to locate those regions which gave 
the largest modifications to the SM results and found, in contrast to the SM, that the 
magnetic amplitude dominates the electric amplitude and that there are large regions of 
the MSSM parameter space for which the magnetic amplitude is greater than that of the 
SM. For regions of high tan/3 (=48) and low (mo,™^) ranging from (225,150) to (275,225) 
the SUSY enhancement of the gluon-mediated exclusive hadronic b decays can be at the 
several percent level and the SUSY penguin processes dominate the SUSY box processes. 
Similar effects of large tan/3 enhancement occur for the photon magnetic dipole operator 
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in b — > S7 and there is significant current interest in its implications for B — > X s 7 and 
B^X s l+l- (l = e,fi,T) 

Calculations of weak decays of hadrons within the framework of the effective Hamiltonian 
involve @j25j: (i) the computation of the quark level decays b —> qq'q' at the electroweak 
scale Mw, (ii) determination of the Wilson coefficients Ct(Mw) through matching of this 
full theory onto the five-quark low energy effective theory in which the W^, t and all SUSY 
particles heavier than My/ are integrated out, (iii) evolution of the Wilson coefficients down 
to the low energy scale fi ~ 0{m^) using renormalization group improved perturbation 
theory |26|| , thereby incorporating the important short distance QCD corrections, and (iv) 
the calculation of hadronic matrix elements for the hadronization of the final-state quarks 
into particular final states, typically evaluated using the factorization assumption |27| . 

In this earlier paper [21J , the QCD corrections arising from renormalization of the short 
distance results down from the electroweak scale to the scale were not included. It was 
argued that these effects were not likely to alter the finding that the magnetic amplitude is 
dominant in the MSSM as the QCD induced mixing effects [2(J produce an enhancement of 
the magnetic dipole operators in the AB = 1 effective Hamiltonian relative to the current- 
current penguin operators associated with the electric monopole amplitude. Furthermore, 
as noted by Gerard and Hou [^], use of the g 2 -dependent SM form factor F[{q 2 ) already 
incorporates the dominant part of the QCD corrections for the current-current penguin 
operators and, therefore, that the main effects of QCD corrections will be the renormalization 
of the strong coupling constant from a s (Mw) to a s (mi,). This would have the effect of 
increasing the penguin decay rates of the SM by the factor rj 2 , where rj = a s (mb) / a s (Mw) ~ 
1.82, and also increasing the MSSM penguin amplitudes relative to those of the SM. 

In the present paper we address each of the stages of the calculation of B decay rates. 
In section II we incorporate the MSSM penguin processes into the effective Hamiltonian 
and fully implement the renormalization-induced QCD corrections. We confirm that the 
MSSM processes are enhanced relative to those of the SM by these QCD corrections and 
that the dominant MSSM processes are indeed the magnetic dipole interactions. In section 
III we compare the quark-level decay rates calculated from the QCD corrected effective 
Hamiltonian with r] 2 T^°\ where r( ) are our weak scale decay rates calculated using q 2 - 
dependent form factors, and find that the Gerard and Hou conjecture is generally valid to 
within 10% for the SM but greatly underestimates the decay rates for the SM+MSSM. The 
hadronization of the quark level processes is studied in section IV where, as an alternative 
to the widely used factorization models, we adopt a simple spectator model. The results 
obtained from this spectator model are given in section V. Section VI presents our discus- 
sions and conclusions. The input numerical values used in our calculations for the CKM 
parameters and quark masses are given in an appendix. 



II. GLUINO PENGUINS AND THE EFFECTIVE HAMILTONIAN 

With the inclusion of gluino-mediated penguin processes from SUSY, the total QCD 
penguin amplitude to lowest order in a s for the decay process b —> q g — ► q q'q' is 

M p eng = _*A [^(pJ^utfaMfyM^vrfa)] (1) 
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where 



% = lAF?tf)PL + F«{f)P R \ + ^[F 2 L (q 2 )P L + F 2 R (q 2 )P R ]. (2) 

Here F\ and F2 are the electric (monopole) and magnetic (dipole) form factors, q = p b — p q 
is the gluon momentum, Pl(r) = (1 T 7s)/2 are the chirality projection operators and 
T a (a = 1, . . . , 8) are the SU(3) C generators normalized to Tr(T a T b ) = \5 ab . In writing (|T]) 
we have neglected the smaller MSSM b — ► qg penguin processes mediated [[F| by charged 
Higgs boson, chargino and neutralino exchanges. The gluino penguin amplitude is enhanced 
relative to that for these processes by both the factor a s /aw and the additional g — g — g 
coupling with its large color factor. 

For the SM, q 2 is not small compared to m 2 u and ml and we must retain the q 2 dependence 



of the u and c contributions to F^(q ). In the limit m q = we then have [pi 

G f 



F^q 2 ) = -|[£ V*V ib h(xi, q 2 ) + V*V tb h{x t )l 

i=u,c 

if (0) = 0, (3) 
F 2 L (0) = 0, 

F 2 R (0) = ^m b £ V*V ib f 2 ( Xl ) (4) 

i=u,c,t 

where V is the CKM matrix and x 2 = rn 2 /My V . Explicit expressions for fi(xi, q 2 ) are given 
in but are not needed in the present study (see later). The q 2 = functions are |]28|j29 



f x ( x ) = — _ [I8x - 29x 2 + 10x 3 + x 4 - (8 - 32x + I8x 2 ) In x] , (5) 

f 2 (x) = 4 ^~ X ^ 4 [2 + 3x-6x 2 + x 3 + 6xlnx]. (6) 

For small x iy / 2 (£j) ~ — \xi whereas f\{xi) ~ — |lnxj. 

For the MSSM form factors we have q 2 <C where mj. is the mass of the loop dj 
squark, and we can use the q 2 = form factors |2l| , [30l , |3lfl 

if (0) = E A^[C a (G)i4(^) + C 2 (it!)£(%)], (7) 
3 

F R (0) = ^ A^^GQA^) + 0,(^)5(5*)], (8) 
x[C 2 (G')C(^)-C 2 (i?) J D(x J )] 

+m 5 A% j [C 2 (G)E(x j ) - 4C 2 (i?)C(^)]}, (9) 

i^(0) = E{KA^ + m 9 A h « Rj \ 
3 

x[C 2 (G)C(%)-C 2 (it:)D(^)] 

+™s A^[C 2 (G)E(^) - 4C 2 (i2)C(%)]}. (10) 
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where m- g is the gluino mass and the coefficient 



describes the rotation from the down-diagonal interaction states to the d mass eigenstates at 
the d — d — g vertices. Here (A,B) are chirality indices, C 2 (G) = 3 and C 2 (i?) = J2 a T a T a = 
4/3 are SU(3) Casimir invariants, j = 1, . . . , 6 labels the d squark mass eigenstates and Xj = 
m| / m|. The matrices V^ L and V^ R are obtained from the (6 x 6) matrix Vj = (Vj L , Vj R ) T 

which diagonalizes the d mass 2 matrix. The MSSM functions are 

1 

6(1 ^x) 4 



A(x) = cn - , 4 [3 - 9x + 9x 2 - 3x 3 + (1 - 3x 2 + 2x 3 ) \nx], (12) 

B(x) = x y l n - lSx + 9x 2 -2x 3 + 6\nx], (13) 

C{x) = ^ [1 - x 2 + 2xlnx], (14) 

= 6 ^~ 1 ^ 4 [2 + 3a:-6x 2 + x 3 + 6xln3:], (15) 

E{x) = — r^[l — x + xlnx]. (16) 

The SUSY masses and diagonalizing matrices needed to evaluate the MSSM form factors 
at M w were obtained (for details see fll9"| ) using two-loop MSSM renormalization group 
equations (RGEs) for the gauge and Yukawa couplings and one-loop MSSM RGEs for the 
other SUSY parameters. The magnitudes of the MSSM form factors satisfy \F^\ > \F±\ > 
|F 2 L | 3> \Ff\ for all regions of the allowed parameter space apart from the narrow region 
tan/? = 2, mi/2 = 150 and niQ > 1000 where \F±\ is slightly smaller than \F% |. Outside 
this region the ratio l-^l/l-^/'l exceeds unity and increases strongly with tan/?. The result 
that Frf is the largest MSSM form factor indicates that, in contrast to the SM, the magnetic 
dipole transition dominates the h decay process in the MSSM. To compare with the SM, we 
note that the ratio of the largest MSSM and SM form factors is |F 2 fl (MSSM)|/|F 1 L (SM)(g 2 = 
0)|<0.4GeV. 

The amplitude (|l|) can be written 

M Peng = M SM + M MSSM ^ 

where 

m sm = -i% v: q v lbh { Xi ^) + V*V tbfl (x t ))0 P 

V2 8 i=uc 

+ *Y E VZVMxfiOs} (18) 



i=u,c,t 



and 



m mssm = ^otsiMwl { h F { {0) o P + F*(0)d P ] + ~[F 2 L (0)O 8 + F 2 R (0)6s}}. (19) 

7T O I 
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Here 



P = A + 6 -^(0 3 + 5 ) (20) 



is a combination of the standard QCD penguin four-fermion operators [26 



3 .5 = (<? a 7 M Uv-A X](^7m^)v t a (21) 

q 1 

4 .6 = {QaYhh-A ^(^9a)v T A, (22) 

9' 

where g' G {w, d, s, c, b}, and 8 is the chromomagnetic dipole operator 

Os = ^m b [q a a^(l + 75)^] G%. (23) 

a and (3 are color indices, the subscripts V=fA represent the chiral projections 1 =F 75 and 
G a is the gluonic field strength tensor. 

In addition to the SM operators, the MSSM gluino penguin processes introduce operators 
of opposite chirality to the O k '- 

o P = d 4 + d 6 -^(d 3 + d 5 ), (24) 



where 



3 .5 = (q a l^b a ) y+A XX^7/«^)v±A (25) 

q' 

4 .6 = (ga7 M Mv+A XX^WJviA, (26) 



and 



8 = £- 2 m b [q a a^(l - l5 )T^] G%. (27) 
The effective Hamiltonian for AB = 1 decays at scale fi = 0{m b ) has the structure 

hm = ^{v^vucMot + cmo%] + v: q v cb [CMoi + c 2 (^)o c 2 } 

-V t *V tb £ [CMO k + CMO k ]} (28) 

fc=3,...,6,8 

where 

0( = (q a q' a h-A % 6/j)v-A, (29) 

o( = (q a q'p)\-A (q'p bah-A, (30) 

(q' G {u,c}), are the tree current-current operators. We omit the photon magnetic dipole 
operator O-j and electroweak penguin operators. The Wilson coefficients C k {n) = Cf M for 
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k = 1, 2 and C k (fi) = Cf M + Cf SSM , C k (fi) = SSM for jfe = 3, . . . , 6, 8 incorporate the 
physics contributions from scales higher than /i and are determined perturbatively at My/ 
by matching to the full SM+MSSM theory. Noting that the u and c contributions from the 
SM penguin to the coefficient of Op in (|i"8"D cancel in this matching process [25], we obtain 
(in the naive dimensional regularization (NDR) scheme) 

Cf«(M w ) = 1-11 (31) 

6 47T 

c su (Mw) = n «^ (32) 

C! M (M W ) = ^1 A k [ fl (x t ) - |], (33) 
Cl M (M w ) = l - f 2 (x t ), (34) 



Cr M (M w ) = (^KV*)- 1 WtiP), 05) 

C 8 MSSM (M W ) = \ (^V; q V tb )-i F 2 L (0) (36) 



where k = 3 6 and 



A fe = (-1,3,-1,3). (37) 



The coefficients {Cf SSM , C'^ SSM } are obtained from {Cf SSM , C 8 MSSM } by the replacement 
^i L 2 (0) - F&(0). 

The Wilson coefficients 

C(M W ) = (d(M w ), . . . , C 6 (M W/ )) T (38) 

evolve under the renormalization group equations (RGE) to 

C(fi) = U 5 ^,M w )C(M w ) (39) 

where U^(ii, Myy) is the five-flavor 6x6 evolution matrix. In next-to- leading-order (NLO), 
U 5 (fi,M w ) is given by [[32}f2B| 



U 5 (fi,M w ) = [1 + ^-J]Ul(fi,M w )[l - ^J] (40) 

47T 47T 

where U®(fi, M w ) is the leading order (LO) evolution matrix and J expresses the NLO 
corrections to the evolution. NLO Wilson coefficients in the NDR renormalization scheme 
were computed from (|39| ) and (f5|) using expressions given in |32|,|26| for J and £7° and the 
two-loop expression for ct s (/i) with five flavors and a 1 f s {Mz) = 0.118 (Aqcd = 0.225 GeV). 

The renormalization group induced mixing between Og and the set {Oi, . . . , Oq} vanishes 
at 1-loop order and it is sufficient to use the LO evolution for Cg(/x): 

CM = rr 14/23 C 8 (M W ) + £ hr ai (41) 

1=1 
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where r] = a s {jj) / a s {Mw) and the constants hi and a\ are tabulated in 

The MSSM generated QCD penguin operators {6 3 (M W ), . . . , 6 6 (M W )} and 8 (M W ) 
renormalize in the same manner as {0 3 (M W ), . . . , 6 (M W )} and 0$(M W ) p4|j . 

The Wilson coefficients Cfc(/x), k = 1, ... ,6,8, for the SM and the combined SM+MSSM, 
together with the Wilson coefficients Ck(f^), k = 3, . . . , 6, for the MSSM are listed in Table 
| for the scales fj, — m& and fi = 2.5 GeV. Inclusion of the MSSM processes generates 
complex Wilson coefficients due to the complex mixing coefficients A b J B j in the MSSM form 

factors. These mixing coefficients arise from the mismatch between the interaction and d 
mass eigenstates at the d — d — g vertices and contain non-removable RGE-generated phases 
in the off-diagonal components of the couplings triggered, in particular, by the complexity of 
the large t-quark Yukawa coupling. The numerical values of J^ABj use d in (^) and (|36D for 



C]f SSM (Mw), k = 3, . . . , 6, 8 are generated from (|TT|) for a given MSSM data set by numerical 
diagonalization of the dj mass 2 matrices at the scale Mw- 
The MSSM results are for the data set 

A = -300 GeV, fj, > 0, tan/5 = 48, m = 275 GeV, m 1/2 = 150 GeV (42) 

which maximizes the effects of SUSY J2]] as measured by the ratio r Pen s(MSSM)/r Pen §(SM) 



of integrated decay rates at scale My/ for b — > qq'q' . This ratio is largest for high tan j3 and 
low (mQ,mi/ 2 ) and for the set ( f4*2| ) is ~ 0.10 for b — > d and ~ 0.085 for b — > s. The 
data set (^) satisfies the constraints imposed |35| on the SUSY parameter space by recent 



experimental bounds on the mass of the lightest chargino (w^± > 91 GeV) and lighter t 
squark (m^ > 80 GeV) and by precision electroweak data constraints on SUSY corrections 
to electroweak parameters but is on the verge of being excluded by the latest measurements 
for Br (B -> X sl ). 

As expected from the observations above on the relative magnitudes of the SM and MSSM 
form factors at scale Mw, the only significant MSSM effects occur in the magnetic dipole 
coefficients. For the data set © we find C 8 MSSM (M W ) = -0.4631 + 0.0224i (b d) and 
—0.4755 — 0.0013i (6 — > s) which are much larger than Cf M (My/) = —0.0953 and comparable 
to the largest SM coefficient Cf M (M w ) = 0.9828. Thus C 8 (fi) is greatly enhanced by the 
MSSM contribution at M w . Although C^ SSM (M W ) is only of comparable magnitude to 
Cf M 6 , it is also enhanced by fl4T| ) and becomes significant. 

The Wilson coefficients in NLO precision depend upon the renormalization scheme 
52],|3B]J3B| and this unphysical dependence is compensated by the perturbative one-loop QCD 



corrections to the matrix elements of the four-quark operators Oi, . . . , O e (and 3 , . . . , 6 ) 
at the scale fi. This is equivalent to using the effective Hamiltonian (|28D with the Wilson 
coefficients Cfc(/i),/c = 1,...,6, replaced by the effective g 2 -dependent Wilson coefficients 
C£ ff (/i, q 2 ) such that 

CM(OM) = Cf^,q 2 )(O k )^. (43) 
The two sets of Wilson coefficients are related by |33| (in the NDR scheme) 

Cf{^ q 2 ) = [1 + ^(r£ + 7 T log ^) kj CM 

+ ^^(Q + c p + c,) (44) 
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where A' k = (0, 0, — 1, 3, — 1, 3) T , the 6x6 constant matrices ry and jv (given in p3| ) 
represent the process- independent parts of the vertex corrections and c t , c p and c g are the 
contributions arising from the penguin-type corrections to the current-current operators 
0\ 2i penguin-type corrections to 3 , . . . , 6 , and the tree level diagram of 8 respectively. 
The quantities c t , c p and c g are given by |J3| 



V+*Vt. 



[V*V ub (- + ^{m u , q 2 )) + V*V cb (- + T x {m c , q 2 ))}, 



tq vtb a 



^iK? 2 )] + (c 4 



3 
C fi 



53 Fi(™<j,<f 

-u,d,s,c,b 



2m b 



where 



and 



T\{m,q 7 



^2 m 2 
= 3 l0g( ^ 



AFA) = -4 / - x) logfl 

Jo 



zx(l 



X 



(45) 
(46) 

(47) 
(48) 
(49) 



Note that Tiijn^q 2 ) is just —fi(x,q 2 ) of pi 



with replaced by [x 2 and is complex for 



q 2 > Am 2 . These strong phases in T\, generated at the uu and cc thresholds, combine with 
the weak CKM phases to produce effective Wilson coefficients C| , k — 3, ... ,6, which differ 
in both magnitude and phase for b and b decays and hence CP violation is manifest. 



Some comments are appropriate about whether or not to include the term c„ in 



This term is absent in most applications of the effective Hamiltonian method prUTS 
was first included by |j33| . The 0$ contribution to (qq'q'\H c s\b) can be written 



.G F Ot s r ,R fr .s'mb r_ 



^)T a u b }[u ql ^T a v q 



and 



(50) 



In the factorization model used by |33|] color considerations preclude q' and q 1 combining 
into the same meson and it is reasonable to assume that p q > = —Pq> in the b rest frame. The 
gluon momentum is then 



Q 



Pb/ m b- 



(51) 



This allows M 8 to be expressed in terms of the penguin operator combination Op given in 
(p0|), and yields the contribution c g to C£ ff , k = 3, . . . ,6. q 2 is conventionally replaced by 
an averaged value (q 2 ) taken to lie in the range ml/A < (q 2 ) < ml/2. If the assumption 



is not valid then c g should be omitted and a proper treatment of the magnetic dipole 
contribution to decay rates must be based upon flSOP. This then involves an integration over 
a range of q 2 determined by kinematical considerations (see later and Refs. [ 21 , Fj ] ) . 

We have computed effective Wilson coefficients with and without the term c g . Although 
the assumption p q > = —p q i is valid for many of the processes (so-called A processes) con- 
sidered in our spectator model (see later), we choose not to make this assumption in our 
effective Wilson coefficients and hence do not include the c g term in our calculations. 



9 



We list in Table [TJ the effective Wilson coefficients Cl s (fi,q 2 ),k = 1,...,6 and 
Cfc ff (//, q 2 ), k = 3, . . . , 6 evaluated from (||) at q 2 = ml/2 and fi = 2.50 GeV for the pro- 
cesses b — > dqq' and b — > dqq' . Table [TII| lists the same quantities for the processes b — > sqq' 
and b — > sqq'. The values of the coefficients are due nearly entirely to the A' k term in 
(fH|). If the input CKM matrix is changed from the Particle Data Group parameterization 
to the Wolfenstein form with the parameter values A = 0.81, p = 0.05 and rj = 0.36 used 
by p5| , and we include c g , we obtain values for the SM C^f which agree very closely with 



those given in |33 



III. DECAY RATES FOR b gg'g 

In calculating the b quark decay rates we allow for the more general case where the 
antiquark produced is not necessarily q'. We firstly consider the SM only and, in accordance 
with most studies of the SM, omit the small magnetic dipole term C| M from the effective 
Hamiltonian ( p8|) in the calculation of decay rates. 

The partial decay rate, b spin averaged and summed over final spin states, has overall 
spherical symmetry. Apart from its overall orientation, a final state is specified by only two 
parameters, say p q = \p q \ and p q > = \p q >\- The partial decay rate in the b rest frame is 

d 2 T G 2 p Pq-Pq' Pq-Pq" , m q ,m q „ 

- -VqPq>E- q „ [ai — — + a 2 —— + a 3 - - ] (52) 



dp q dp q i IT 3 E q E q/ E q E q » E q >E q ,. 



where 



and 



a x = \dx + d 2 + d 3 + d^ 1 + 2\d x + d 4 \ 2 + 2\d 2 + d 3 \ 2 , 
a 2 = |rf 5 + c / 6 | 2 + 2|d 5 | 2 + 2|rf 6 | 2 , 

a 3 = 3fc[(3di + d 2 + d 3 + 3d 4 )d* 6 + (d x + 3d 2 + 3d 3 + d 4 )d* 5 ] (53) 



<f) = Vq>> q Vq> b CJCa*. f), 4,..., 6 (/i, g 2 ) = ~V*V tb C^Jp, q 2 ) . (54) 



Rerepg-pgi = [m 2 b + m 2 q ,, -m 2 -m 2 , -2m b E ql/ ]/2, p q -p q » = [m 2 b + m 2 q , -m 2 -m 2 q „ -2m b E q >]/2 
and the angles between the particle velocities must be physical, for example | cos(# g(? ')| < 1. 
Note that, for q e {d, s}, only the tree current-current terms d\ )2 contribute when q' ^ q" e 
{u, c} whereas, for AC = AU = transitions, only the penguin terms d 3t ^$ contribute for 
q' = q" G {d, s} and both tree and penguin terms contribute for q 1 = q" G {u, c}. 
For massless final state quarks, the decay rate is |37[ 



r = i2r f de[ ttl (e) + a 2 (fl][i - f + & ( 55 ) 

Jo o 2 3 

where r = G 2 F ml / 192n 3 and ^ = q 2 jm\. For SM pure penguin decays, the QCD corrected 



decay rates at scale /i calculated from ( 55|) generally agree to better than 30% (see Table [Tv| 



with r] 2 T^ where T^ ^ are the weak scale decay rates calculated using g 2 -dependent form 



factors [21] and the factor rj accounts for the running of a s from M\y to \i. To this extent 



the results support the conjecture of Gerard and Hou [24] mentioned in the introduction. 
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To study the effects of the MSSM on the quark decay rates, we note that the only 
significant effects of the MSSM are confined to the magnetic dipole coefficients Cg(/i) and 
C 8 (fj,) (see Table |) and the terms Ck, k = 3, . . . , 6 of the effective Hamiltonian ( |28| ) can 
be neglected in any calculation of decay rates. Since \C$(fi)\/\C 8 ([i)\ < 0.2 for both b — > d 
and b — > s transitions, it would be reasonable to also neglect the Cg term so that, if the 
assumption (|5lD about the momenta of q' and q" is applicable, the major effects of the 
MSSM could be incorporated in a modified value of c g . However, for the general case, the 
c g term should be omitted from C| , k = 3, . . . , 6 and the magnetic dipole terms directly 



included. For massless final state quarks the dipole contributions add the terms ETJ 



l(\ds\ 2 + N 2 )(^ - \ + f ) - 1 + d 4 + - O 2 (56) 



to the integrand of (p>oD, where 

dM = -V*V tb C 8 (/i), dM = -V*V tb CM, (57) 

thus requiring the imposition of a lower cutoff on q 2 (taken to be 1.0 GeV 2 ). The MSSM 
gluino penguin processes so included are greatly enhanced relative to the SM penguin am- 
plitudes by the renormalization-induced QCD corrections (see Table |V|) and the resulting 
decay rates F for the SM+MSSM far exceed the ^ 2 -scaled weak scale decay rates T . We also 
note from Table |V| that, whereas the QCD corrections only mildly affects the magnitudes 
of the CP asymmetries acp = (T — T)/(T + T) for the SM, they can decrease the magnitude 
and change the sign of acp for the SM+MSSM. 

The impact on b — > sg of new physics exhibited by an enhanced chromomagnetic coupling 



Cg involving an unconstrained new CP phase has recently been studied by Hou [14| within 
an effective Hamiltonian framework. g 2 -dependent decay rates and CP asymmetries are 
given for phase differences of tt/4, tt/2 and 37r/4 between the Cg amplitude and the complex 
SM C 3v .. i6 penguin amplitudes. In all cases Hou finds large rate asymmetries above the 
cc threshold, a feature not present in our present calculations, or in our earlier weak scale 
calculations incorporating a g 2 -dependent SM F±(q 2 ) form factor [21], because of the much 
smaller MSSM phases (for b — > s the phases of the dominant MSSM form factors and 



Fi are nearly constant [^TJ at ~ —0.016 over the allowed MSSM parameter space). 



Although the QCD corrected SM+MSSM decay rates are quite sensitive to the choice 
of lower cutoff for q 2 , it is clear that the values obtained for the MSSM data set ( ^2|) are 
unacceptably large. For example, assuming a maximum branching ratio of 1% for b — > sdd 
gives an experimental upper bound on the b — > s decay rate of ^ 3 x 10~ 15 GeV, a factor 
of more than 10 smaller than our calculated rate. Even though the MSSM data set used in 
the present calculations is not typical of the MSSM parameter space, having been chosen to 
maximize the hadronic b decay rate at My/, our findings here suggest that experimental data 
on hadronic b decays will exclude a similar low (^0,^1/2) region of the MSSM parameter 
space to that excluded by the b — > 57 constraints f35fl . 

From here on we focus on decay rates in the SM for b decays which are not pure penguin 
processes. We therefore neglect the magnetic dipole interaction and retain final state quark 
masses. If we make the common assumption that the coefficients a^, i = 1,2,3 are constant, 
evaluated at q 2 = ml/2 for example, then the decay rate is 
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r = r M£ = i)/! + « 2 (e = \)h + « 3 (e = |)/ 3 ] (58) 

where the phase space integrals Jj,z = 1,2,3 have to be computed numerically. Decay 
rates computed from fl58|) are shown in Table M, together with branching ratios for various 
b decays. Branching ratios are less sensitive to the value chosen for rrib and have been 
obtained by dividing the decay rates by the total decay rate r Total computed using the 
0(a s ) QCD-corrected semileptonic decay rate 

2(T 

F(6 -> d-*) = V,\V cb \ 2 h [1 - -^/(p c )], (59) 

D7T 



where p c = (m c /mb) , with f(p c ) = 2.51 |38| . 



IV. A SPECTATOR QUARK MODEL FOR B -> h x h 2 

We now wish to make the transition from decay rates at the quark level to decay rates for 
two body hadronic decays B — > hih 2 - Calculation of the required hadronic matrix elements 
(hih 2 \Ok\B) from first principles is currently not possible and several approximate schemes 
are available in the literature. Much work has been done on the factorization model ]27| and 



the more recent generalised factorization model |H|^,(|l|] of this process in which final state 



interactions are neglected and the hadronic matrix elements are factorized into a product of 
two hadronic currents: {hi\J\ ^\B) (h 2 \J2 |0). The operators 0\ and O^q are Fierz transformed 
into a combination of color singlet-singlet and octet-octet terms and the octet-octet terms 
then discarded. The singlet-singlet current matrix elements are then expressed in terms 
of known decay rates and form factors. Consequently, the hadronic matrix elements are 
expressed in terms of the combinations 

«2fc-l = Cfk-l + Jj-Czki a 2fc = C*2fc + -^-Qfc-i (60) 

where k = 1, 2, 3 and the number of colors N c is usually treated as a free phenomenological 
parameter in order to compensate for the discarded octet-octet terms. Factorization works 
reasonably well for B decays to heavy hadrons in that a\ and a 2 seem to assume 



universal values, and it is argued |33[] to also account for B decays to light hadrons such as 



B — > Kn and B — > -kit. The generalised factorization approach has been criticised p5|j25 



because the effective Wilson coefficients fl44| ) are gauge dependent, implying that the value 
of N c extracted from comparisons of factorization model predictions with data cannot have 
any physical meaningQ These authors present an alternative model based upon Wick 
contractions in the matrix elements of the NLO effective Hamiltonian. Furthermore, the 
applicability of factorization to the hadronization of outgoing light quark pairs moving in 



opposite directions has been questioned ||7|] and it has been emphasised there that hadronic 



B decays, especially as they are originated by three partons, are phase space driven. 



1 Gauge invariant and infrared finite effective Wilson coefficients have recently become available 



m 



12 



In the spectator model |42| alternative to factorization, the b quark and spectator quark 
are treated as quasifree particles with a distribution of momenta due to their Fermi motion 
relative to the B meson. In these models the b quark is treated as a virtual particle with 
an invariant mass mb(p s ) satisfying m 2 {p s ) = m 2 B + m 2 — 2w,bJp1 + m s- We, however, 



consider here a tentative spectator model based upon the idea of duality between quark and 
hadron physics at the high energies of b quark and B meson decays. The decays at the quark 
level, even including the penguins, are basically short distance processes. In our proposed 
spectator quark model the long distance hadronization is largely a matter of incoherently 
assigning regions of the final quark phase space to the different mesonic systems. 

For example, we consider a B meson bu or bd to be a heavy stationary b quark accompa- 
nied by a light spectator constituent antiquark which has a spherically symmetric normalized 
momentum distribution P(\p s \ 2 )d 3 p s . The total meson decay rate through a particular mode 
is then assumed to be 

T =f -T-J- P(\Ps\ 2 )d 3 p s dp q dp q ,, (61) 

J dpgdpq' 

equal to the initiating decay rate. Neglecting, for the moment, any constraints due to quark 
color, we suppose the spectator antiquark to combine with the quark i(i = q' or q) to form 
the meson system. For example, if i = q' we assign a mass M q i s to the system such that 

M q's = iPq' + Ps) ■ {Pq' + Ps) = m\ + m 2 + 2{E q ,E s - p q ,p s COs9 q > s ). (62) 

Constraining p q > and p s to have mass M q / S , we can infer from ( j6T|) that 

dT r d 2 T 

= 2M q , s J — — - P(\p s \ 2 )5(M 2 , s - m\, -mr 8 - 2(E q ,E s - p q ,p s cos 9 q > s ))d 3 p s dp q dp ql . 

(63) 



Hence 



where the integration region is restricted by the condition | cos^^l < 1. 

We also assign a mass M qq u to the second quark antiquark system such that 

M qq" = (Pq + Pq") ' (Pfl + Pf) = (P& ~ Pq') ' (Pb ~ Pq') = m b + m \' - 2m h E q ,. (65) 

The variable E q i, and hence p q ', determines the mass M qq n. Taking this mass to be the 
independent variable, we have 

d 2 T 2nM q , s M qq -„ f E q ,p s d 2 T 2 



dM q , s dM qqll m b J p 2 , dp q dp, 

We call this mode of quark and antiquark combination process A. Finally, by integration, we 
compute the partial decay rates T(M q i s , M qq >r) into quark systems with masses less than M q > s 
and M qq ". With suitable binning we equate these partial decay rates with corresponding 
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rates into mesons. To calculate an exclusive decay into a particular two meson system would 
require one to follow the flow of quark spins. Although the model could be generalized to 
do so, we attempt to take large enough mass bins to enclose all the spin combinations. 
When searching for CP asymmetries we would advocate looking at more inclusive decay 
rates which will be larger and hence more readily measurable. Of particular interest are the 
quark antiquark systems forming the lowest mass 0~ and 1~ meson states as data exists for 
charmed quark systems which can be used to test the spectator model. Also, data which 
should exhibit matter-antimatter asymmetry is eagerly awaited for the rare decays into light 
quark systems. 

It is also possible that the spectator antiquark combines with the quark q, for which we 

get 



dM qs dM q ,qu m b J p\ dp q dp, 

We call this process B. 

In some meson decays, for example B~ — ► (ud) + (uu) which is initiated by the quark 
decay b — > udu, the spectator u could have combined with the u or the d. However, for 
light meson systems such as 7r~p°, we find that the different combinations come from very 
different regions of phase space of the initiating b decay and we conclude that we are not 
double counting. 

Turning now to the flow of color, we examine what may be regarded as two extreme 
possibilities. In the first, called here model I, we take the formula (0), (|66|) and flBTD at face 
value, that is we make no attempt to follow the flow of color and assume that all color flow 
is taken care of by the gluon fields in the meson system. In the second, called here model II, 
we consider the possibility that the lowest mass meson states are only formed if the quark- 
antiquark pair are in a color singlet state. That is, we assume that the color disruption 
caused by a quark-antiquark pair in a color octet state will result in more complex meson 
systems than the lowest mass 0~ or 1~ states. 

Projecting out the color singlet states results only in a modification of the coefficients 
q;i 2 of For process A represented by (|66|) they become 



'3 

whereas for process B represented by flBTD they are 



ai = 3\di + irf 2 + ^rd 3 + d A \ 2 

«2 = 3|-d 5 + d 6 \ 2 , (68) 



1 i 

ctx = 3|-di + d 2 + d 3 + -cLi| 2 

a 2 = 3|rf 5 + ^ 6 | 2 - (69) 
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V. SPECTATOR MODEL RESULTS 



Calculations have been performed for two spectator quark momentum distributions. One 
is a fixed spherical bag model |48| in which the spectator is confined in a spherical cavity of 



radius R and has a wave function of the form 



= 1 /AsinfAO 
7T V 2 r 



This implies a momentum distribution 



n\P,?) = (-) 3 , tl'7'( A L 

7T |p s | 2 (A 2 - |Ps| 2 ) 2 



and a mean square momentum |p s | 2 = A 2 . The parameter A is determined from the mean 
square B meson radius 

^=t4-^ < 72 > 

Taking Vr 2 = 0.55 fm yields the value A = 0.6 GeV which we have used in our calculations. 
As a check on the sensitivity of our results to the model, we have also used a Gaussian 
spectator quark distribution with \p s \ 2 = A 



2. 



P(|Ps|2) = ( 27A )3/2e • (73) 

For the maximum mass of the quark- ant iquark systems we take a value midway between 
the lowest mass 1~ state and the next most massive meson. Thus we take, for (cu) or (cd), 
M ca = 2.214 GeV between the D* (2.007) and the D* (2.420); for (cs), M cS = 2.323 GeV 
between ££(2.112) and £> s *(2.533); for (su) or (sd), M sa = 1.081 GeV between 7^(0.892) 
and (1.270); for (ud), M ud = 0.877 GeV between p(0.770) and a (0.984); and for (uu) and 
(dd), M U u = M d g = 0.870 GeV between cj(0.782) and 7/(0.958). For the (cc) system we make 
an exception and take the maximum M c5 = 3.552 GeV between x(3.417) and if)(2s) (3.686). 



The results are presented in Table [VT] and compared, where data is available, with the 
sum of the branching ratios into mesons with masses less than the above cutoff masses. All 
calculations are based upon partial decay rates (|52|) with the a, coefficients evaluated at 
q 2 = ml/2. The semileptonic decay rates are also shown in Table |VT[ For these decays the 
lepton momenta cover all the allowed phase space. 

We have repeated the calculations using the SM+MSSM effective Wilson coefficients 
of Tables |TJ and [III] but neglecting the very small C eS coefficients. The MSSM effects are 
insignificant for the b — > udu transitions. For type A processes all branching ratios are 
slightly increased (~ 1%). For b — > usu transitions the MSSM effects are more significant 
and result in branching ratio decreases for type A decays of 23% for b — > usu and 28% for 
b — > usu. The b — > dsd and b dsd branching ratios are also decreased by about 30% 
Inclusion into our spectator model of the MSSM-enhanced magnetic dipole term is a matter 
for future investigation. 
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VI. DISCUSSION AND CONCLUSIONS 



We have investigated two aspects of hadronic B decays. The first is the effects of 
renormalization-induced QCD corrections on the inclusion of SUSY gluino penguin con- 
tributions to the hadronic b decays b — > qg — > qq'q', q' G {d,s}, which have been studied 
by calculating the Wilson coefficients of the effective SM+MSSM Hamiltonian. The SUSY 
enhancement of these gluon mediated processes at M\y by the MSSM magnetic dipole tran- 
sition, which is especially significant for [^TJ large tan/3 and low (mg,™^), is further am- 
plified by the QCD corrections, resulting in quark decay rates at scales \i ~ 0{m^) for pure 
SM+MSSM penguin processes which greatly exceed those obtained by simply scaling the 
rates at Mw by the factor rj 2 to allow for the renormalization of the strong coupling con- 
stant from Mw to fi. The rates are also significantly larger than estimated experimental 
bounds, suggesting that experimental data on hadronic b decays will exclude a similar low 
(mo, mi/2) region of the MSSM parameter space to that excluded by the b —* s-y constraints 
|[35|| . We find that, whereas the QCD corrections only mildly affect the magnitudes of the 
CP asymmetries acp, they can change their sign. 

The second aspect of hadronic B decays studied is the hadronization of the final state 
quarks. We have adopted a spectator model as an alternative to the widely used factorization 
models. The semileptonic decay rates given in Table |VI| indicate semi-quantitative agreement 
between our model and data. Where there is data, the same is true for the hadronic decays 
that proceed by process A, suggesting that color suppression does not have a large influence 
on these processes. This is not so for process B hadronic decays which we find to be very 
sensitive to color suppression and there is clear evidence for color suppression in the data. In 
fact, for decay modes which can proceed only through process B, the data are for the most 
part upper bounds. The decays B° -> ^(ls) + (K°, K* °) and B + -> if>(ls) + (K+,K* + ), 
for which the experimental branching ratios are (2.24 ±0.3) x 10~ 3 and (2.46 ±0.37) x 10~ 3 
respectively, are process B decays in the context of our model. Taking a cutoff mass M c5 = 
3.257 GeV to include ip(ls) but exclude the x(l-P) we find, with no color suppression, a 
branching ratio of 9.22 x 10~ 3 . This suggests a color suppression factor of about four, not 
the orders of magnitude obtained from insisting that the cc pair are produced as a color 
singlet (model II). 

The decays into light quark systems exhibit matter-antimatter asymmetry. The present 
calculations indicate that, averaged over a few low mass 7r mesonic states, the B + and B° 
branching ratios are larger than those for B~ and B°; the opposite is true if a K meson is 
involved. The calculations suggest the asymmetry is not large. In the effective Hamiltonian 
used here, this asymmetry comes only from the "strong" phases of the penguin diagrams. 
The renormalization group construction of the effective Hamiltonian, although incorporating 
important QCD improvements, does not include the "strong" phases which must be present 
even at the first order of a s corrections. The penguin includes just one 0(a s ) correction to 
the "strong" phase; it would be of interest to include all 0(a s ) corrections to the "strong" 
phase. 

We have compared our results for the light quark systems with the factorization model 
calculations of |36j which presents results for pseudoscalar and vector meson final states. 
Because of the constraints of spin and isospin, which are averaged over in our model, we 



are reluctant to use cuts which bracket individual particles. However, the results of [36 
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when summed over the pseudoscalar and vector final particles are generally smaller than 
our calculations, but by less than factors of two. Such factors can easily be accomodated 
within the uncertainties of just the CKM parameters for the light quark couplings. 
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APPENDIX: INPUT PARAMETERS FOR NUMERICAL CALCULATIONS 

We use the standard Particle Data Group j|3] parameterisation of the CKM matrix with 
the central values 

012 = 0.221, 0i3 = 0.0035, 23 = 0.041 (Al) 

and choose the CKM phase S13 to be tt/2. 

Following Ali and Greub |33|], we treat internal quark masses in penguin loops as con- 
stituent masses with the values 

m d = mu = 0.2 GeV, m s = 0.5 GeV, m c = 1.5 GeV, m b = 4.88 GeV. (A2) 

For the light quarks in the spectator model we take 

m u = 0.005 GeV, m d = 0.01 GeV, m s = 0.2 GeV (A3) 

and a spectator quark mass of 0.01 GeV. 
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TABLES 



TABLE I. Wilson coefficients C k (fi) for the SM and the combined SM+MSSM and Wilson 



coefficients 


Ck for the MSSM at the renormalization scales \x - 


= m b and fi = 2.50 GeV. The MSSM 


results are 


for the MSSM data set 


©• 






olvl 






SM+MSSM 






Z.OU \j& V 




2.50 GeV 


C\ 


1 .U 1 O 1 


1.1ZDD 


1. 0767+0. OOOOi 


1.1266+0. OOOOi 


Co 


U.loll 


— U.Z 1 01 


-0.1811+0.0000i 


-0.2751+0. OOOOi 


c 3 


0.0119 


0.0181 


0.0119+0.0000i 


0.0182+0. OOOOi 


c 4 


0.0331 


-0.0461 


-0.0332+0.0000i 


-0.0462+0. OOOOi 


c 5 


0.0094 


0.0113 


0.0094+0.0000i 


0.0113+0.0000i 




0.0398 


-0.0598 


-0.0399+0.0000i 


-0.0599+0.0000i 


c 8 


0.1449 


-0.1585 


-0.4768+0.0160i 


-0.4556+0. 0144i 


cf 






-0.39(-8)+0.19(-9)i 


-0.43(-8)+0.21(-9)i 


Cf 






0.71(-8)-0.34(-9)i 


0.68(-8)-0.33(-9)i 


Cf 






-0.14(-8)+0.65(-10)i 


-0.79(-9)+0.39(-10)i 


Cf 






0.11(-7)-0.55(-9)i 


0.13(-7)-0.65(-9)i 


C s 






-0.0768+0.0000i 


-0.0976+0.0000i 


a Numbers 


in parentheses denote 


powers of 10 





TABLE II. Effective Wilson coefficients Cf(fi,q 2 ) for the SM and the combined SM+MSSM 
and effective Wilson coefficients Cf for the MSSM at the renormalization scale \i = 2.50 GeV and 
momentum q 2 = m 2 /2 for b —* d and b — ► d transitions. Results are given for the MSSM data set 

(ID- 







SM 


SM+MSSM 




b^d 


b^d 


b^d 


b^d 


cf 


1.1679+0.0000i 


1.1679+0. OOOOi 


1.1679+0.0000i 


1.1679+0. OOOOi 


cf 


-0.3525+0.0000i 


-0.3525+0.0000i 


-0.3525+0.0000i 


-0.3525+0.0000i 


cf 


0.0233+0.0018i 


0.0250+0.0047i 


0.0188+0.0019i 


0.0204+0.0046i 


cf 


-0.0540-0.0053i 


-0.0588-0.0142i 


-0.0403-0.0057i 


-0.0451-0.0137i 


cf 


0.0172+0.0018i 


0.0189+0.0047i 


0.0127+0.0019i 


0.0143+0.0046i 




-0.0675-0.0053i 


-0.0723-0. 0142i 


-0.0539-0.0057i 


-0.0587-0.0137i 


cf 






-0.0010+0.0000i 


-0.0010+0.0000i 


cf 






0.0029+0.0000i 


0.0029+0.0000i 


cf 






-0.0010+0.0000i 


-0.0010+0.0000i 


cf 






0.0029+0.0000i 


0.0029+0.0000i 
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TABLE III. Effective Wilson coefficients Cf(n,q 2 ) for the SM and the combined SM+MSSM 
and effective Wilson coefficients Cf for the MSSM at the renormalization scale \i = 2.50 GeV and 
momentum q 2 = m 2 /2 for b — > s and b — > s transitions. Results are given for the MSSM data set 

ft 



SM SM+MSSM 
b — > s b — > s b — > s b — > s 



Cf 


1.1679+0.0000i 


1. 1679+0. OOOOi 


1.1679+0.0000i 


1.1679+0.0000i 


cf 


-0.3525+0.0000i 


-0.3525+0.0000i 


-0.3525+0.0000i 


-0.3525+0.0000i 


r<e& 


0.0249+0.0032i 


0.0248+0. 0030i 


0.0202+0.0032i 


0.0201+0.0030i 


Cf 


-0.0585-0.0095i 


-0.0583-0.0090i 


-0.0445-0.0095i 


-0.0442-0.0090i 


Cf 


0.0188+0.0032i 


0.0187+0.0030i 


0.0141+0.0032i 


0.0140+0.0030i 


cf 


-0.0721-0.0095i 


-0.0718-0.0090i 


-0.0581-0.0095i 


-0.0578-0.0090i 


cf 






-0.0010+0.0000i 


-0.0010+0.0000i 


cf 






0.0031+0.0000i 


0.0031+0.0000i 


cf 






-0.0010+0.0000i 


-0.0010+0.0000i 








0.0031+0.0000i 


0.0031+0.0000i 



TABLE IV. Comparison of QCD corrected integrated decay rates T(/j,) at scale fj, calculated 
from (|55|) and (^) for pure penguin b — > q and b —> q, q G {d,s} transitions with the scaled 
decay rates T^°\fj,) = rj 2 T^ (M\y) where T^°\M\y) are the weak scale decay rates calculated using 
g 2 -dependent form factors and rj 2 = (a s (fi) / 'a s (M\y)) 2 accounts for the renormalized a s . Results 
are calculated for the cutoff q^ n = 1.0 GeV 2 and for the MSSM data set (|42"|). The numbers in 
parentheses denote powers of 10. 



r(GeV) b -» 


d b^ 


SM 
d b^ 


s 


b^ 


s b^ 


SM+MSSM 
d b^d b^ 


s b- 


TW(M w ) 1.480(- 


-17) 2.210(- 


-17) 


2.991(- 


-16) 


2.921(- 


-16) 2.062(- 


-17) 2.926(- 


-17) 4.080(- 


-16) 3.997( 


r(°)(m 6 ) 4.898(- 
T(m b ) 6.304(- 


-17) 7.311(- 
-17) 9.076(- 


-17) 
-17) 


9.897(- 
1.303(- 


-16) 
-15) 


9.667(- 
1.276(- 


-16) 6.823(- 
-15) 2.931(- 


-17) 9.681(- 
-15) 2.825(- 


-17) 1.350(- 
-15) 5.054(- 


-15) 1.323( 
-14) 5.065( 


r(°)(2.50) 7.594(- 
r(2.50) 8.103(- 


-17) 1.134(- 
-17) 1.222(- 


-16) 
-16) 


1.534(- 
1.714(- 


-15) 
-15) 


1.499(- 
1.674(- 


-15) 1.058(- 
-15) 2.665(- 


-16) 1.501(- 
-15) 2.542(- 


-16) 2.093(- 
-15) 4.558(- 


-15) 2.051( 
-14) 4.570( 
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TABLE V. SM phase space integrals Jj, i = 1,2,3, decay rates and branching ratios B for 
various b decays calculated using effective Wilson coefficients C| ff (//, q 2 ) evaluated at the renormal- 



ization scale fj, = 


2.50 GeV. The hadronic decay rates T are 


calculated from (|58f) 


using q 2 fixed at 


m 2 /2. The numb 


ers in parentheses denote powers of 10. 








Decay process 


h h h 




T (GeV) 


B 


b — > c d u 


0.5011 




1.846(-13) 


5.007(-l) 


b — ► c s u 


0.4911 




9.280(-15) 


2.517(-2) 


b — > c d c 


0.1765 0.1766 8.69(- 


-2) 


2.921(-15) 


7.706(-3) 


b — > c s c 


0.1702 0.1706 8.49 - 


-2) 


5.507(-14) 


1.494(-1) 


b —* u d u 


1.0000 0.9952 8.38(- 


-6) 


2.354(-15) 


6.386(-3) 


b — > u d u 






2.586(-15) 


7.014(-3) 


b — > u s u 


0.9884 0.9823 8.34(- 


-6) 


2.462(-15) 


6.677(-3) 


b — > iis u 






2.228(-15) 


6.043(-3) 


b — > tide 


0.4982 




6 362 C— 17) 


1.726(— 4) 


b — > tt d c 






6.362(-17) 


1.726(-4) 


6 — ► u s c 


0.4882 




1.235(-15) 


3.351(-3) 


6 — > c e~ v e 


0.5012 




4.991(-14) 


1.354(-1) 


b^ c n~ 


0.4982 




4.961(-14) 


1.346(-1) 


b — > C T~ v T 


0.1121 




1.117(— 14) 


3.030(-2) 
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TABLE VI. Experimental and spectator model branching ratios for B meson decays. 
The notation, for example, B° — > (D~,D*~) + (tt + ,p + ) indicates the branching ratio 



B° -> [D~ (1.869) or D* "(2.010)] + (vr+ or p+). The 


process classification A or 


B refers to 


whether the spectator antiquark combines with the quark q', Eqn. 


(|66|) or q, Eqn. 


(67). The 


model classification I or II refers to no color suppression and total color suppression respectively 


in the formation of the decay mesons (see the text). For the decays into light quark systems, the 


data are for the most part upper bounds which are not shown. The mc 


del predictions are in accord 


with these bounds. Only b quark decays are shown, as in 


the PDG tables. 




Decay Experiment 


Process 


Model I 


Model II 


b — > c l + v\ 








B + -> (D°,D*°) + l + vi (7.16 ± 1.3)xl0~ 2 




8.6xl0~ 2 




B° -» (D~,D* -) +l + Vi (6.6±0.5)xl0~ 2 




8.6xl0~ 2 




7 , l 

b — > u r v\ 








B° -» (ir-,p-) (4.3 ± 1.6)xl0~ 4 




5.2xl0~ 4 




b — > c u d 








B° -» (D-,D* -) + (vr+, p+) (2.04 ± 0.53) x 10~ 2 


A 


3.00xl0~ 2 


2.73xl0~ 2 


B+ -» (D°, D* °) + (7T+, p+) (3.88 ± 0.58) x 10~ 2 


A 


3.00xl0~ 2 


2.78xl0~ 2 




B 


1.39xl0~ 2 


1.55xl0~ 5 


£° (D°,D*°) + (7r ,77,p ,u;) < 3.2xl0~ 3 
b — > c s c 


B 


1.39xl0~ 2 


1.55xl0~ 5 


5+ -> (Z)°, 5* °) + (D+, L>* +) (6.1 ± 2.3) x 10~ 2 


A 


3.2xl0~ 2 


2.9xl0~ 2 


B° — > (£>-,£>* -) + (£>+,£>*+) (4.8 ± 1.8) xlO" 2 


A 


3.2xl0~ 2 


2.9xl0~ 2 


B° -► (V-(ls), X(l-P)) + 0^°, ^* °) < 6x 10~ 3 


B 


3.7xl0~ 2 


4.1xl0~ 5 


B+^^(ls), X (lP)) + (K + ,K*+) 


B 


3.7xl0~ 2 


4.1xl0~ 5 


b ^ u du 








B ^ (ir+,p+) + (n-,p-) 


A 


1.07xl0~ 4 


0.95xl0~ 4 


B+^(ir+ t p+) + {*°,ri,p ,u,) 


A 


1.07xl0~ 4 


0.95xl0~ 4 




B 


1.09xl0~ 4 


1.9xl0~ 7 


B° — ► (7T°, r/, w) + (vr°, r/, w) 


B 


1.09xl0~ 4 


1.9xl0~ 7 


6 — ► it d u 








(7r+,p+) + (7r-,p-) 


A 


0.95 xl0~ 4 


0.87xl0~ 4 


B~ -» (^-,p-) + (^ ,r ? , / 9°, W ) 


A 


0.95 xl0~ 4 


0.87xl0~ 4 




B 


0.97xl0~ 4 


1.9xl0~ 7 


5°^ (^,pV) + (7r ,r,,pV) 


B 


0.97xl0~ 4 


1.9xl0~ 7 


6 — ► u s -u 








£?° — > (-RT+, K* + ) + p~) 


A 


0.503xl0~ 4 


0.422 xl0~ 4 


B+^(K+,K*+) + (<jr ,r,,p ,u;) 


A 


0.503xl0~ 4 


0.422 xl0~ 4 




B 


1.09xl0~ 4 


0.93xl0~ 6 


B° ^ (K°,K* ) + (-ir , V ,p°,Lj) 


B 


1.09xl0~ 4 


0.93xl0~ 6 


b — > u s u 








B° ^ (K-,K*-) + (7r + ,p+) 


A 


0.624 xl0~ 4 


0.531 xl0~ 4 


B- -> (^-,X*-) + (vr ,r / ,p°, W ) 


A 


0.624 xl0~ 4 


0.531 xl0~ 4 




B 


1.21xl0~ 4 


0.93xl0~ 6 


£?° — » (^T°, °) + (7T°, r/, p°, u) 


B 


1.21xl0~ 4 


0.93xl0~ 6 



d s d 
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B+ - (K ,K* ) + (7r+,p+) 
6 — > <i s J 

5-^(^ ,K* ) + (7T-,p-) 

5° -► (K°,K*°) + (^ ,r/,p ,^) 



0.491 xlO" 4 0.410xlO~ 4 

0.491 xlO" 4 0.410xlO~ 4 

1.082xl0 -4 0.920xl0~ 6 

0.499xl0~ 4 0.417xl0~ 4 

0.499xl0" 4 0.417xl0~ 4 

1.096xl0 -4 0.914xl0~ 6 
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